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Abstract 

The temporal evolution of combustion fiowtields estab- 
lished by the interaction between wedge-shaped bodies and 
explosive hydrogen-oxygen-nitrogen mixtures accelerated to 
hypersonic speeds in an expansion tube is investigated. The 
analysis is carried out using a fully implicit, time-accurate, 
computational fluid dynamics code that we developed 
recently for solving the Navier-Stokes equations for a chem- 
ically reacting gas mixture. The numerical results are com- 
pared with experimental data from the Stanford University 
expansion tube for two different gas mixtures at Mach num- 
bers of 4.2 and 5.2. The experimental work showed that flow 
unstart occurred for the Mach 4.2 cases. These results are 
reproduced by our numerical simulations and, more signifi- 
cantly, the causes for unstart are explained. For the Mach 5.2 
mixtures, the experiments and numerical simulations both 
produced stable combustion. However, the computations 
indicate that in one case the experimental data were obtained 
during the transient phase of the flow'; that is, before steady 
state had been attained. 

Introduction 

The development of new detonation-wave-based propul- 
sion devices, such as the ram accelerator [ I ]. oblique detona- 
tion wave engine [2] and pulse detonation engine [3], has 
renewed interest in studying the combustion of premixed 
gases flowing at hypersonic speeds. Recently, a new experi- 
mental technique — based on the expansion tube — was devel- 
oped by Srulijes et al. [4] and Morris et al. [5,6] to 
investigate such flows. In the expansion tube, combustible 
gas mixtures are accelerated to hypersonic speeds and then 
allowed to interact with stationary test bodies. The resulting 
flowfield can produce various modes of shock-induced com- 
bustion, ranging from decoupled shock-deflagration waves 
to overdriven detonation waves. 

The expansion tube consists essentially of a single tube 
divided into three or four sections by diaphragms, as shown 
in Fig. 1, which is a schematic of the classic expansion tube 
developed by Trimpi [7]. When the test gas comprises an 
explosive mixture, a buffer zone (not shown in Fig. 1) con- 
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taining an inert gas is added between the test gas and the 
driver gas, to prevent autoignition following rupture of the 
primary diaphragm [4]. 

The principal advantage of the expansion tube over the 
shock tube and shock tunnel in applications involving pre- 
mixed combustible gases is that the operating cycle does not 
involve stagnation of the test gas. Therefore, explosive mix- 
tures can be accelerated to hypersonic velocities without 
autoignition. 

Other methods have also been used to study 
shock-induced combustion flowfields. For example. Lehr [ 8] 
used ballistic range experiments, in which projectiles were 
fired at high speeds into premixed hydrogen-air mixtures. 
However, one disadvantage of this technique is that measure- 
ments are difficult to make on a moving projectile. In con- 
trast, in the expansion tube the test models are stationary, 
thereby simplifying data collection. 

The main difficulties associated with the expansion tube 
(and pulse facilities in general) are the short test times avail- 
able (up to half a millisecond) and the relatively low' total 
pressures attainable. A rule of thumb w hen air is the test gas 
is that, to assure flow' establishment, the pulse facility must 
supply a “slug'’ of air that is three times longer than the test 
model [9,10]. For chemically reacting flows, especially if 
recirculation zones and shock wave/boundary layer or deto- 
nation wave/boundary layer interactions are present, the slug 
of test gas may have to be longer than the recommended 
value, as shown in this paper. Thus, because of the finite test 
time available, it may not be possible to determine experi- 
mentally the required flow' establishment time for a given 
reactive mixture and operating condition. 

Computational fluid dynamics (CFD) simulations can 
help resolve the above issue, because they do not have test 
time limitations. Thus important questions, such as the exist- 
ence of a steady state and the test time required for its estab- 
lishment, can be answered by CFD. Also. CFD can be used 
to extrapolate test results to the high pressures where propul- 
sion systems will actually operate and which cannot cur- 
rently be attained in the laboratory. Thus CFD can serve as a 
powerful method for complementing experimental work. 
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In this paper we present results of our numerical simula- 
tions of five experiments conducted in the Stanford Univer- 
sity expansion tube by Morris et a l [5,6] and make 
qualitative comparisons with their observations. Recently, 
Choi et al.[ 11] also carried out numerical simulations of two 
of the problems examined in this work. Although we do not 
show their results, we make comparative observations. The 
main objective of our study is to investigate the details of 
combustion initiation and temporal evolution of the flow- 
field, in order to better understand the flow physics and 
explain the experimental observations. In addition, these 
experiments serve as new benchmark test problems for our 
ongoing program on developing and validating an accurate, 
efficient and robust CFD code for studying high speed, 
chemically reacting, viscous flows [ 12-16]. 

Expansion Tube Cycle 

A brief explanation of the operation of the expansion 
tube is given below, in order to put into perspective the time 
domain of our calculations. Figure 1 shows a schematic dia- 
gram of the expansion tube cycle. Following rupture of the 
primary diaphragm at time t = /q, a primary shock wave (.V| 
in Fig. 1 ) propagates into the test gas, and an expansion wave 
into the driver gas. On reaching the end of the driven section, 
the primary shock ruptures the secondary diaphragm and a 
secondary shock wave (a^) propagates downstream into the 
expansion section, while an expansion wave moves upstream 
into the test gas. This expansion wave is washed down- 
stream, since the gas in region 2 is moving at supersonic 
speeds. Test time begins with the arrival of the test gas/accel- 
erating gas contact discontinuity (C 2 ) at the model, and ends 
with the arrival of the expansion w ave. The state of the gas in 
region 5 determines the test condition. 


Q = J 'lP|. P 2 > -P/v P". P*’. ef (2) 

The dependent variables are the mass density of the ith 
species p, (/ = 1 ,...,/V), the velocity components u and v and 
the total energy per unit volume e. J is the grid Jacobian and 
F and G are the inviscid flux vectors in the £ and rj direc- 
tions, respectively. F v and G v are analogous viscous fluxes. S 
and S v are axisymmetric source terms and W is the chemical 
source term. A detailed description of the terms in equation 
(1 ) and additional state and constitutive equations needed for 
system closure are given by Yungster [17]. 

All five cases examined in this study involve mixtures of 
hydrogen, oxygen and nitrogen. The chemical reaction 
mechanism was taken from Jachimowski [18], except reac- 
tions involving nitrogen are not considered, because they 
become important only at higher Mach numbers than consid- 
ered here. The resulting combustion mechanism consists of 
19 elementary reactions among 8 reacting species and the 
inert species No [12]. 

Numerical Method 

The numerical method used for solving equation (1) is 
described by Yungster and Radhakrishnan [12] and summa- 
rized below. For clarity in presentation, only the 
two-dimensional Euler equations are considered 
here — extension to the viscous case is, however, straightfor- 
ward [12]. The equation set is discretized using the back- 
ward differentiation formula (BDF) method [19]. because its 
accuracy and efficiency in solving the ordinary differential 
equations arising in combustion chemistry have been well 
established (20,21 1. 


Our interest in the present work is to examine the (low 
and combustion processes resulting from the interaction of 
the high speed combustible mixture with the test model. 
Therefore, our simulations are performed for time / > tj (see 
Fig. 1 ). 


Define the time step At* 1 = ^ +1 - A the time step ratio 
r r = Aft/Af 1 ' 1 and the change in the vector of dependent vari- 
ables AQ n = Q ,l+] - Q n . Then the second-order, variable-step 
BDF method applied to equation (I) gives for the time inter- 
val [t, r t fl+{ ] 


Numerical Formulation 

Governing Equations 


AQ 


./< * 


yA Q"~ k ' -VAi"lF J+W 2. k -Fj_ l/2 , k + 
G j. k + 1 /2 - Gj k 1 /2 - W j' j. 1 


(3) 


The conservation form of the unsteady Navier-Stokes 
and species transport equations describing two-dimensional 
or axisymmetric, chemically reacting flows can be written in 
general curvilinear coordinates (£, rj) as follows: 


30 3(F-F V ) 3(G-G V ) 
3/ + 3^ + 3n 


+ j(S-s r ) = w 


( 1 ) 


w here t is time, the parameter j is zero for two-dimensional 
flow and one for axisymmetric flow and Q is the vector of 
dependent variables for a gas mixture containing N species: 


where the BDF method coefficients y and (3 are 


y = 


1 + 2 r t 


p = 


1 + r, 

1 + 2r, 


(4) 


The terms F and G are the numerical fluxes in the £ and 
rj directions. They are computed using Yee’s second order 
total variation diminishing (TVD) scheme [22]. Equation (3) 
is then linearized in a conservative manner and solved itera- 
tively, by using a lower-upper relaxation procedure consist- 
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ing of successive Gauss-Seidel (LU-SGS) sweeps. At each 
time step the iterative process of producing successively 
improved approximate solutions to equation (3) is continued 
until a suitable convergence criterion is satisfied 1 12,19]. 

The inversion of large matrices is avoided by partitioning 
the system into reacting and nonreacting parts. Conse- 
quently, the matrices that have to be inverted are of the same 
size (NxN) as those that arise in the commonly used point 
implicit methods. A fully coupled interaction between the 
flow and chemistry is, however, maintained through the sub- 
iteration process. An important advantage of the present 
method is that, because it is fully implicit, it is stable for 
large values of the CFL number, thereby enabling the use of 
relatively large time steps to minimize computational cost. 

Results 

In this section we present results of our numerical simu- 
lations of reaction initiation and development in the Stanford 
University expansion tube [5,6] for five different cases. We 
also compare our results with available experimental OH 
planar laser induced fluorescence (PLIF) data and, although 
not shown here, with the numerical solutions of Choi et al 
[11], who recently conducted a CFD study of the expansion 
tube flowfield. The test condition and reactive gas mixture 
composition are given in Table 1 for all five cases. In this 
table p , T and M refer to the pressure, temperature, and Mach 
number, respectively. The subscripts 5, 10 and 20 denote the 
expansion tube regions shown in Fig. 1 and the subscript CJ 
refers to the Chapman-Jouguet condition corresponding to 
the gas mixture at state 5; A/ C j was computed with the 
NASA chemical equilibrium code CEA [23]. Note that all 
five cases involve flow at superdetonative speeds; that is, the 
test gas velocity is greater than the CJ detonation speed. 

Although helium was used as the accelerating gas in the 
experiments, hydrogen was selected in the numerical simula- 
tions for computational simplicity. This change in the accel- 
erating gas does not significantly affect the results in the time 
interval of interest, because the two gases have similar spe- 
cific heat ratios and the accelerating gas is not involved in the 
combustion process. This simplifying assumption was also 
made by Choi et al [ II]. 

Figures 2a-c show the experimental OH PLIF images for 
cases 1-3. Cases 4 and 5 resulted in flow' unstart, and no OH 
PLIF image was given by Morris et al. [5,6]. 

Case 1 

The first case studied the flow at Mach 5.2 (i.e., — 5.2) 

of a 2H 2 + 0 2 + 1 7N 2 mixture over a 40° wedge-shaped pro- 
jectile of length L = 30 mm. as shown in Fig. 3a. For this test 
problem we assumed turbulent flow and adiabatic projectile 


wall. Zero-order extrapolation was used for the top and right 
boundaries. The Baldwin-Lomax turbulence model [24] was 
utilized, and a 313x85 grid was found to adequately resolve 
the flow features. 

Figures 3a-h show the temporal evolution of the reacting 
flowfield. by means of water mass fraction contours (top half 
of each figure) and normalized temperature (77Tjo) contours 
(bottom half of each figure). These figures show the reflec- 
tion of the expansion tube's secondary shock (sy see Fig. 1) 
by the wedge (Fig. 3a). followed by arrival of the contact dis- 
continuity (cy see Fig. 1 ) and high-speed combustible gas 
mixture (Fig. 3b). As the combustible mixture reaches the 
projectile, an oblique shock wave begins to form (Fig. 3c). 
Subsequently, at t ~ 29 (is, combustion is initiated along the 
boundary layer (Fig. 3d). Combustion is initially limited to 
the boundary layer, but it then begins to expand laterally at a 
location just upstream of the projectile shoulder (Figs. 3e-h). 
At t ~ 92 (is the solution becomes essentially time indepen- 
dent (Fig. 3h); that is, steady state is obtained. 

The lateral expansion of the combustion front (Figs. 
3e-h) can be explained as follows. Combustion in the bound- 
ary layer causes the oblique shock wave to curve slightly, 
thereby increasing its strength. Consequently, the induction 
time decreases along the length of the wedge and combus- 
tion intensifies near the projectile shoulder. However, the 
shock wave-combustion front interaction is weak, and so a 
decoupled shock-deflagration wave system is obtained. 

The computed steady state solution is also shown in 
Fig. 4, but in the form of OH mass fraction contours. Com- 
paring this figure to the corresponding experimental OH 
PLIF image (Fig. 2a) shows the good qualitative agreement 
between the computed structure of the combustion zone and 
the experimental result. At t- 92 jls, the approximate flow 
establishment time for case 1 , the test “slug” of combustible 
gas mixture is 166 mm long, or about 5.5 times wedge 
length — significantly longer than the suggested value 
(~ 3xL) for nonreacting test gases [9,10]. 

Cases 2-5 

Cases 2-5 consider hypersonic flow over two 20° wedges 
arranged symmetrically, as shown in Fig. 5a. This configura- 
tion is similar to the hollow projectiles used by Thibault et 
al. [25] and Sasoh et al. [26] in their ram accelerator studies. 
The test flow domain is 61 mm long, with a constant 
cross-sectional area region 12.7 mm high. After Choi et a\. 
[11], who carried out numerical simulations of cases 2 and 3 
(M 5 = 5.2; see Table 1 ), we assume adiabatic walls, laminar 
flow and zero-order extrapolation for the outflow boundary, 
in order to facilitate comparison with their results. For these 
two cases also, the flow features were resolved adequately 
with a 313x85 grid. 
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The first stage of flow development, that is. from arrival 
of the secondary shock until combustion initiation, is very 
similar for cases 2-5. Therefore, we show this period of flow 
development (which covers approximately the first 75-80 (is) 
only for case 2. which involves the flow' of a 2 H 2 + O 2 + 
17N2 mixture at Mach 5.2 (see Table 1). The results are 
given in Figs. 5a-l, in the form of normalized temperature 
(7>T|o) contour lines. In Fig. 5 and subsequent figures (6-9). 
the contour range, given in the form (minimum, maximum, 
increment), applies for the plot in which the range is speci- 
fied and for every subsequent plot, until a new r range is 
defined. 

Time is measured from the instant when the secondary 
shock (,v 2 ) is 3 mm upstream of the leading edge of the test 
model (Fig. 5a). As the secondary shock moves over the test 
model, it is reflected and a complex system of oblique 
shocks and expansion waves is established (Figs. 5b-f). After 
approximately 32 ps, the oblique shocks coalesce into a nor- 
mal shock that moves slow ly upstream (Figs. 5g and 5h). At 
approximately 40 (is, just before the normal shock wave 
moves out of the computational domain (Fig. 5h), the contact 
discontinuity (r 2 ) (followed by the high-speed test gas) over- 
takes the normal shock (Fig. 5i). (It should be pointed out 
that in the Stanford University experiments [5.6] the reacting 
mixture arrives at the projectile at a later time than we have 
assumed. We chose an earlier time, in order to avoid having 
to use an unnecessarily long computational domain. This 
simplifying assumption should not significantly affect the 
subsequent reacting flow. The same approach was adopted 
by Choi et al \ 11).) 

The interaction between the contact discontinuity and the 
normal shock produces a new r normal shock (Fig. 5i) that 
travels downstream, and can be seen leaving the computa- 
tional domain at t ~ 56 (is (Fig. 5j). As the high speed test 
gas mixture begins to flow over the model, a shock wave sys- 
tem begins to form that is attached to the two wedges (Fig. 
5 j ). This system consists of a pair of weak and strong oblique 
shocks interconnected by a normal shock. At approximately 
70 (is the two strong oblique waves collide and the normal 
shock disappears (Fig. 5k). Later, at t ~ 76 jus. this interac- 
tion results in the formation of a new r normal shock between 
the two weak oblique shocks waves at the leading edges 
(Fig. 51). At this time the first stage of flow development 
(i.e., the preignition regime following secondary shock 
arrival at the projectile) is essentially over, because combus- 
tion initiation occurs soon, as discussed next. 

The temperature behind the normal shock is high enough 
to cause ignition of the mixture. Combustion begins at t ~ 85 
(is, as shown in Fig. 6a, which gives normalized pressure 
(pip | 0 ) contour lines overlaid on OH mass fraction contours. 
Production of OH is observed in Fig. 6a to occur down- 
stream of the normal shock, separated by a short induction 
zone. Subsequently, the flame front moves closer to the nor- 


mal shock (Figs 6b and 6c). and concurrently the width of 
the normal shock decreases: eventually, the oblique-normal 
shock system transitions into a regular oblique shock reflec- 
tion (Fig. 6d). The temperature behind the reflected oblique 
wave is, however, not high enough to maintain combustion, 
and so the flame detaches from the oblique shock (Figs. 6e 
and 60- The flame is then washed downstream, out of the 
computational domain (Figs. 6g-i). (Transient combustion 
phenomena in expansion tubes is not surprising, inasmuch as 
it was observed in our previous studies of unsteady reacting 
flowfields in these devices [ 1 3].) 

Combustion in the boundary layer is observed to begin at 
t ~ 103 Jls, mainly around the region w'here the oblique 
shock impinges on the body surface (Fig. 6d). The shock 
wave-boundary layer interaction creates a separated flow 
region, which acts as a flameholder (Figs. 6d-i). When 
steady state is attained at t ~ 197 (is, combustion persists 
only inside the boundary layer upstream of the projectile 
shoulder, within the separated flow region and along the reat- 
tached boundary layer downstream (Fig. 6i). This time corre- 
sponds to an effective mixture slug 319 mm long, or 5.2 
times model length — again much longer than the recom- 
mended value for air [9,10]. 

The structure of the computed steady state flowfield (Fig. 
6i) shows good qualitative agreement with the experimental 
OH PLIF image (Fig. 2b). Note in particular that the experi- 
mental and numerical results both show that combustion is 
restricted to the separated flow' region and boundary layer 
downstream — reaction is not observed in the central flow 
region. (It should be pointed out that the area imaged by the 
PLIF system did not extend to the model shoulder (see Ref. 6 
and Fig. 2b). Hence some of the reacting region in Fig. 6i is 
not visible in Fig. 2b.) 

We also compared our results with those reported 
recently by Choi et al [1 1] (results not shown here). Excel- 
lent agreement between the two solutions was obtained, both 
at steady state (compare Fig. 2c in [ 1 1 ] and Fig. 6i) and dur- 
ing flow establishment (i.e,, the transient phase; see Fig 6 in 
[11] and Figs. 6a-i). 

Case 3 also involves flow at Mach 5.2. but a different 
mixture (2H 2 + 0 2 + 12N 2 ) than in case 2 (see Table 1 ). The 
temporal evolution of the computed postignition flowfield is 
shown in Figs. 7a-i, as superimposed OH mass fraction con- 
tours and normalized pressure (p/p\o) contour lines. The 
flow development during this stage is similar to that of case 
2: an initial period of combustion in the core region, fol- 
lowed by flame blow' out, and then establishment of a stable 
reaction region in the separated flow region and in the 
boundary layer upstream and downstream of this region. But 
in case 3 combustion in the core flow' region persists for a 
longer time period than in case 2 (compare Figs. 6h and 7f, 
which were both generated at approximately the same time 
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point), due to the higher energy content of test gas mixture 3 
(because of the smaller N 2 concentration in mixture 3 than in 
mixture 2: see Table 1 ). 

The flame remains attached to the oblique shock until 
approximately 132 ps (Fig. 7f). Then detachment of the 
flame occurs and the combustion zone moves 
downstream — eventually out of the computational domain 
(Figs. 7g-i). The steady state flowfleld (Fig. 7i). attained after 
approximately 194 (is, is similar to that for the previous case 
(see Fig. 6i), but the separated flow region is somewhat 
larger. The effective gas mixture slug length required for 
reacting flow establishment (i.e., the length corresponding to 
t - 194 ps) is 318 mm, or again approximately 5.2 times 
model length. 

Comparing the calculated steady state flowfleld structure 
for case 3 (Fig. 7i) with the corresponding experimental OH 
PLIF image (Fig. 2c) shows that the experimental result indi- 
cates combustion along the center of the flow domain, 
whereas the numerical solution does not. However, the com- 
puted solution at earlier times, before steady state establish- 
ment, does display core flow reaction, as described above 
(see Figs. 7b-7h). Moreover, careful examination of Figs. 
7e-g shows that the experimental image given in Fig. 2c cor- 
responds to the numerical solution obtained at some time 
between -120 ps and -140 ps (but not at steady state. Fig. 
7i). In order to ascertain the reason for the discrepancy at 
steady state between the computed and experimental flow- 
fields, we compared our results with those of Choi et al[ 11). 
As in case 2. excellent agreement was noted between the two 
computed flowflelds, both at steady state (not shown here; 
see Fig. 2d in [ 1 1 ]) and during flow development (see Fig. 5 
in [11]). Also. Choi et al. [II) remark on the absence of 
combustion in the central region of the computed flowfleld, 
in contrast to the experimental result. It can therefore be con- 
cluded that the experimental image shown in Fig. 2c was 
obtained during the transient phase; that is, before steady 
state had been attained. In other words, the test time was not 
long enough for the reacting flow to become fully estab- 
lished 

Cases 4 and 5 consider the same test model and gas mix- 
tures as cases 2 and 3, respectively, but the test gas Mach 
number is 4.2 (see Table 1 ). For both cases 4 and 5, Morris et 
al. [6] reported flow unstart, sometimes followed by a deto- 
nation wave traveling up the expansion tube. No OH PLIF 
image was given for either case, and the cause of flow unstart 
could not be determined conclusively. 

The postignition results of our numerical simulations for 
cases 4 and 5 are given in Figs. 8 and 9, respectively, as over- 
laid OH mass fraction contours and normalized pressure 
( p/p | o ) contour lines. Our initial calculations indicated that 
the 313x85 grid used for cases 1-3 was not fine enough to 
capture accurately the unstart process. To estimate the 
required axial grid spacing, we calculated the length of the 


induction zone behind a Mach 4.2 normal shock wave, by 
using the NASA kinetics code LSENS [21,27). The induc- 
tion zone lengths were 1.1 mm for case 4 and 0.8 mm for 
case 5, respectively. A new 508x180 grid was created, with a 
variable grid spacing of 0.18 mm at the leading edge and 
0.05 mm at the shoulder. Therefore, the induction zone 
would be resolved with a minimum of four grid points. In 
order to attain the higher grid resolution, the test model 
length was reduced to 42 mm. 

Figures 8a-j show the reacting flow' development for case 
4 and ;>75ps. As in case 2, the intersection of oblique 
shocks results in the formation of a normal shock (see Figs. 
6a and 8a). However, note that the normal shock is now 
located farther upstream than in case 2. Combustion in the 
core region is observed to exist from t~ 75 (is to / ~ 127.5 
(is, when the flame detaches from the reflected oblique shock 
and is washed downstream, leaving the computational 
domain (Figs. 8a-c). At t ~ 160 (is, there is no longer any 
combustion in the central flow' region (Fig. 8d). 

Boundary layer combustion begins at approximately 
100 (Lis (Fig. 8b), and is clearly observed at 127.5 (is (Fig. 
8c). A separated flow' region with intense combustion is 
formed, and the separated flow- produces new weak oblique 
shocks (Fig. 8d). This region continues to grow', and propa- 
gates upstream of the projectile shoulder, generating addi- 
tional weak oblique shocks (Figs. 8e-j). At t - 190 [is. the 
multiple oblique shocks produced by the separated (and 
reacting) flow increase the temperature in the core region to 
above the autoignition temperature, thereby giving rise to a 
new combustion front in this region (Fig. 8f). 

An oblique detonation wave forms, which propagates 
upstream of the projectile shoulder (Figs. 8g-h). thus increas- 
ing the pressure in that region and accelerating upstream 
flame propagation in the boundary layer. At t - 219 (is, 
boundary layer combustion reaches the projectile leading 
edge (Fig. 8h). At this time, the boundary layer combustion 
cannot propagate further upstream; instead, it begins to prop- 
agate laterally towards the center of the flow channel (Fig. 
8i). Subsequently, the laterally propagating flame interacts 
with the upstream-moving oblique detonation wave, forming 
a normal detonation wave, which is apparent at t ~ 262 (is 
(Fig. 8j). Note that in this case (and the next), because of det- 
onation, the peak pressures are much higher than in cases 2 
and 3 (see Figs. 6-9). 

At this time, the calculation was stopped, because we had 
captured the essential features of the experimental observa- 
tions — flow unstart and formation of a normal detonation 
wave that propagates upstream of the test model. 

Figs. 9a-j illustrate the results of our numerical simula- 
tion of case 5, which involves a more energetic mixture than, 
but the same test gas pressure and temperature as. case 4. 
Figs. 8a-j and 9a-j show that the reacting flowfleld develop- 
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ment is very similar for the two cases, but in case 5 combus- 
tion in the core region is more intense than in case 4 and this 
flame does not leave the computational domain, However, as 
in case 4, the oblique shocks generated by the separated flow 
reintensify combustion in the central region (Fig. 9f). Note 
that at this time (/ ~ 196 ps) combustion in the boundary 
layer has not propagated as far upstream as in case 4 (com- 
pare Figs. 8f and 9f). 

Again, as in case 4, an oblique detonation wave is formed 
that propagates upstream of the projectile shoulder (Figs. 
9g-h), with an attendant rise in pressure, and flow unstart 
ensues. Also observed is the upstream propagation of the 
boundary layer flame, which reaches the leading edge of the 
projectile (see Figs. 9h and 9i), and then expands laterally 
into the core region (Fig. 9i). In this case, the oblique detona- 
tion wave propagating upstream splits into three distinct 
cells (Figs. 9h and 9i). The interaction between the oblique 
detonation wave system and the laterally expanding reacting 
boundary layer transforms this complex wave system into a 
normal detonation wave (Fig. 9j), as in case 4 (see Fig. 8 j). 

Conclusions 

Numerical simulations of combustion flowhelds resulting 
from the interaction of reactive mixtures accelerated to 
hypersonic speeds with projectiles in an expansion tube 
facility were performed, and results presented for five differ- 
ent cases. In all cases our computations reproduced qualita- 
tively the experimental observations. 

In the first case, which involved hypersonic flow around a 

40° wedge, stable shock-induced combustion, wherein the 
flame front was decoupled from the oblique shock, was 
obtained. The computed steady state combustion structure 
was in very good agreement with the experimental OH PLIF 
image. Combustion was shown to originate in the boundary 
layer. A slug of test gas mixture 5.5 times projectile length 
was required to attain steady state; that is, for reacting flow 
establishment. This length is about 80% longer than the rec- 
ommended value for air. Thus additional research is required 
to determine flow' establishment time for reacting flows in 
pulse facilities. 

For the symmetrical 20° wedge configuration and test gas 
Mach number of 5.2 (cases 2 and 3), the numerical solution 
reached steady state at approximately 190 jlis. This time cor- 
respond to a test gas length of 5.2 times projectile length, 
again significantly longer than the recommended length for 
air. 

In cases 2 and 3, the computed steady state showed com- 
bustion only in the boundary layer, especially within the sep- 
arated flow' region created by the shock wave/boundary layer 
interaction. We found that our results agreed w'ell with the 
numerical solutions of Choi et al[\\]. For case 2, the experi- 


ments showed the same combustion structure as the calcula- 
tions. However, for case 3, the experiments showed 
combustion along the central region of the flow, while the 
numerical results did not. This difference was attributed to 
the experimental data being recorded during the transient 
phase of the reacting flow; that is. before steady state had 
been established. 

For the second configuration and test gas Mach number 
of 4.2 (cases 4 and 5), the simulations predicted flow unstart 
and the generation of a detonation wave, in agreement with 
the experimental observations. The flow unstart was shown 
to be caused by the boundary layer combustion propagating 
upstream of the projectile shoulder and towards its leading 
edge. This flame propagation produces a series of oblique 
shock waves that reignite the core flow', creating an oblique 
detonation wave whose interaction with the laterally expand- 
ing boundary layer flame gives rise to a normal detonation 
wave that propagates upstream. 
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Table 1: Test condition and reactive mixture composition for cases 1-5. 
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Accelerating Gas (H 2 ) 

Test Gas 
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0.018 
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0.1 13 
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2H 2 + 0 2 + I7N 2 
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Figure I . Schematic diagram of expansion tube cycle. Numbers identify flow regions as defined by Trimpi [7], 
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Figure 4. Computed steady state OH mass fraction contours Figure 3. Temporal evolution of expansion tube flowfield 

for case 1 (from Yungster and Radhakrishnan [15]). f° r case 1* The top half of each figure gives water mass frac- 

tion contours and the bottom half normalized temperature 
(77r 10 ) contours (from [15]). 
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(d) t — 1 8.42 (is (e) f = 21.89 ns (f) / = 26.58 p.s 



(j) / = 55.80 [is (k) f = 69.61 |is;7'/7' 10 (1.0, 5.9,0.18) (1) / = 76. 18 (J.s 


Figure 5. Normalized temperature (TIT i0 ) contour lines showing development of initial (i.e., preignition) flowfield for case 2. 
(In Fig. 5 and subsequent figures (6-9), the contour range, given in the form (minimum, maximum, increment), applies for the 
plot in which the range is specified and for every subsequent plot, until a new range is defined.) 


NAS A/CR— 1999-209304 


10 






p/p 10 (2, 381,24) 



N AS A/CR— 1999-209304 


12 


Figure 9. Temporal evolution of reacting fiowfield for case 5, showing overlaid 
OH mass fraction contours and normalized pressure (pip\§) contour lines. 
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